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ABSTRACT 

A simple physical model for long-duration gamma ray bursts (GRBs) is used to fit the redshift (z) 
and the jet opening-angle distributions measured with earlier GRB missions and with Swift. The effect 
of different sensitivities for GRB triggering is sufficient to explain the difference in the z distributions 
of the pie-Swift and Swift samples, with mean redshifts of (z) = 1.5 and (z) = 2.7, respectively. 
Assuming that the emission properties of GRBs do not change with time, we find that the data can 
only be fitted if the comoving rate-density of GRB sources exhibits positive evolution to z > 3 - 5. 
The mean intrinsic beaming factor of GRBs is found to range from « 34 - 42, with the Swift average 
opening half-angle (8j) ~ 10°, compared to the pie-Swift average of (8j) ~ 7°. Within the uniform 
jet model, the GRB luminosity function is oc L~ 3 - 25 , as inferred from our best fit to the opening 
angle distribution. Because of the unlikely detection of several GRBs with z < 0.25, our analysis 
indicates that low redshift GRBs represent a different population of GRBs than those detected at 
higher redshifts. Neglecting possible metallicity effects on GRB host galaxies, we find that 1 GRB 
occurs every 600,000 yrs in a local L* spiral galaxy like the Milky Way. The fraction of high-redshift 
GRBs is estimated at 8 - 12% and 2.5 - 6% at z > 5 and z > 7, respectively, assuming continued 
positive evolution of the GRB rate density to high redshifts. 
Subject headings: gamma-rays: bursts — cosmology: theory 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are brief flashes of 7 rays 
occurring at an average rate of a few per day throughout 
the universe. The ultimate energy source of a GRB is be- 
lieved to be associated with a c atastrophic event assoc i- 
ated with black-hole formation (jMeszaros fc Reeall997D . 
This event takes place through the collapse of the core 
of a massive star in the case of long-duration GRBs, and 
due to merger or accretion-induced collapse events for the 
short-hard class of GRBs. For the class of long-duration 
GRBs, which is the only type considered in this paper, 
liberation of energy from the GRB source is thought to 
result in a very high temperature, baryon-dilute fireball 
that expands to reach highly relativistic speeds and form 
a collimated outflow. Dissipation of the directed kinetic 
energy of the relativistic blast wave through internal and 
external shocks is the standard mod el for producing the 
prompt radiation and afterglow fsee lMeszarosll2006l for 
a recent review). 

Evi dence of jetted GRBs is obtai ned from ra- 
dio dWaxman, Kulkarni fc Fraill [1998) and opti- 
cal (jStanek et all 119991) observations of achromatic 
breaks in the afterglow light curves. The structure 
of these jets is, however, still an open question. The 
two leading models are (1) the uniform jet model, in 
which the energy per solid angle d£*/d£l is roughly 
constant within a well-defined jet opening angle 9j, 
but drops sharply ou tside of 9j, with 0, d if fering 
between GRBs (see IMeszaros. Rees fc Wiiersl 119981: 
IMeszaros fc Reeslll997l): and (2) the univer s al structured 
jet model (jRossi. Lazzati fc Reesl 120021 : IZhang et al.l 
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120041 ) . in which all GRB jets are intrinsically identical 
with the directional energy release d£*/d£l dropping 
approximately as the inverse square of the opening angle 
fr om the jet a xis. 

iFrail et all (|2001h and iPanaitescu fc Kumar! (|2001l ) 
found that the absolute 7-ray energy releases, after cor- 
recting for the jet opening angles inferred from the af- 
tergl ow light curv e s, are clustered near £* 7 = 5 x 10 50 
ergs. iBloom et al.l (|2003[ ). treating a larger sample, con- 
firmed this clustering around £* 7 ~ 1.3 X 10 51 ergs (the 
total energy release £* = § dQ (d£*/dQ) > £* 7 depends 
on the assumed effici ency for 7-ray energy con version). 
In a later analysis by iFriedman fc Blooml (|2005f ) . still in 
the pie-Swift era, a broader spread by about one order 
of magnitude is found in the beaming corrected 7-ray 
energy release around a mean value of £» 7 w 2 x 10 51 
ergs (see Fig. 1). Since we use jet opening angles de- 
rived from these analyses, which are based on the uni- 
form jet model, we also adopt this mod e l in o ur treat- 
ment. The sample of lFriedman fc Blo om ( 2005}) is taken 
to provide a good representation of the pie-Swift ob- 
servations of GRBs. GRBs can in princip le exist (e.g., 
ISchmidtl [19991 IMeszaros fc Meszarosl 119961) and be de- 
tected to very high redshift, z > 10 (jCiardi fc Loebl l200d: 
lLamb fc Reichartl feOOO) . which holds promise that GRBs 
can be used to probe the early universe. Their utility 
for this purpose depends on the frequency and bright- 
ness of high-redshift GRBs, which in turn depends on 
the rate-density evolution of GRB sources and the evolu- 
tion of GRB source properties with cosmic time. Given 
the assumption that the GRB rate density is propor- 
tional to the measure d star format i on rate (SFR) his- 
tory of the universe (|Totanil [1991 E ain fc Nataraianl 
20001) . the redshift distribution of GRBs can be predicted 
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Fig. 1 . — Histogram of the corrected absolute 7-ra y energy release 
£, 7 for pre-Swift GRBs l|Friedman fc Bloonj|2005l) . 

(|Bromm fc Loebl [20021 ) . Correlations between spectral 
properties and energy releases may also permit GRBs 
to be used to determine cosmological parameters (e.g., 
IGhirlanda et alj|2004t iMortsell fc Sollermanll2005f ). 
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Fig. 2. — The cumulative redshift distribution of GRBs 
for 41 GRBs in the pre-Swift sample (thin histograms; 
I Friedm an & Bloom 2005), an d 16 G RBs in the Swift sample (thick 
histograms; Jakobsson ct al. 20061). The median redshifts of the 
pie-Swift and Swift bursts are z ta 1.25 and z pa 2.70, respectively. 

Wi th the launch of the Swift satellite (|Gehrels et al.l 
2004), rapid follow-up studies of GRBs triggered by 
the Burst Alert Telescope (BAT) on Swift became 
possible. A fainter and more distant population of 
GRBs than found with the pie-Swift satellites CGRO- 
BATSE BeppoSAX INTE GRAL, and HETE-2 is de- 
tected (jBerger et al.l l2005t ). The mean redshift of 
41 pre-Swift GRBs that also have measured beam- 
ing breaks (|Friedman fc Blooml l2005f) is (z) — 1.5, 
while 16 GRBs discovere d by Swift have (z) = 
2.72 |Jakobsson et all l200ll (see Fig. 0). iBagolv et all 
(2006) performed five different statistical tests to 
compare the redshift distributions of Swift and pre- 
Swift GRBs assuming, as the null hypothesis, that they 
are identical distributions. They conclude that the red- 
shift distributions are different, that is, all five tests re- 
ject the null hypothesis at the > 97% significance level. 
Thus, the pre-Swift and Swift data sample distinct sub- 
sets of the overall GRB population due to different de- 
tector characteristics and followup capabilities. 

Here, we develop a physical model to understand the 
differences between the redshift distributions of Swift and 



pre-Swift GRBs, taking into account the different detec- 
tor triggering properties. The GRB model is based on 
the uniform jet model for GRBs, though the typical prop- 
erties are assigned in accordance with GRB observations. 
The GRB rate density and distribution of jet opening an- 
gles are the principle unknowns, with parameter values 
for the mean intrinsic duration and absolute 7-ray en- 
ergy release of GRBs determined from observations and 
best fits to the data. A GRB is detected if its flux ex- 
ceeds the flux threshold within the energy window of a 
detector ( the actual detecti on criteria are more compli- 
cated; see Band 2003, 2006, and discussion below). For 
a given rate density of GRBs, the redshift and opening 
angle distributions are calculated and fit to the data. 

We consider only data for the long, soft population 
of GRBs. The separation of the long duration GRBs 
from the short, hard GRBs, which clearly comprise 
a separate source population is complicated both by 
the extended emission observed with Swift from GRBs 
that would normally be cla ssified as short hard GRBs 
(e.g.. lBarthelmv et al.l l2005). and the appearance of the 
unusual GRB 0606 1 4 that defies simple classification 
(IGehrels et al.1 [20061: [Gal- Yam et all [20061 : lFvnbo"et~aTI 
2006). We therefore additionally assume that the data 
set used to test our model represents only, or mostly, 
classical long duration GRBs. We return to this point in 
the discussion. 

In § 2, we develop the formulation for bursting sources, 
and we fit the data in §3. Inferences for the differences 
between the Swift and pre-Swift redshift distributions 
are also made in § 3, including implications for beaming 
factors and low-z GRBs. We conclude in § 4 by discussing 
the derived GRB source rate density, the nature of the 
low-redshift GRBs with weak energies, and the predicted 
fraction of high-redhift GRBs. Derivation of the bursting 
rate of sources in a flat A cold dark matter cosmology is 
given in the Appendix. 



2. STATISTICS OF COSMOLOGICAL GRBS 

2.1. Model GRB 

We approximate the spectral and temporal profiles of 
a GRB occurring at redshift z by an emission spectrum 
that is constant for observing angles 9 = arccos /i < 6j — 
arccos/ij to the jet axis during the period Ai»; thus the 
duration of the GRB in the observer's frame is simply 
At = (1 + z)At* (stars refer to the stationary frame, and 
terms without stars refer to observer quantities). The 
vF v spectrum is therefore written as 

vF v ee fXt) £ f tpk S{x)H{n; N ,l) H(t;0,At) , (1) 

where the spectral function S(x) = 1 at x = e/e p k — 
e*/tpk*- At observer time t, e = hv/m e c 2 is the di- 
mensionless energy of a photon in units of the electron 
rest-mass energy, and e p k is the photon energy at which 
the energy flux f e takes its maximum value f tpk . The 
quantity i/(/i;/ij, 1) is the Heaviside function such that 
H(fi; fJtj, 1) = 1 when fj,j < fi < 1, that is, the angle 
of the observer with respect to the jet axis is within the 
opening angle of the jet, and H(p; fij, 1) = otherwise. 
One possible approximation to the GRB spectrum is a 
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broken power law, so that 



S(x) 



\ e pkJ 



H{e pk - e) + 



\£pk 



H( 



£pk) 



(2) 

where the Heaviside function H (u) of a single index is de- 
fined such that H(u) = 1 when u > and H (u) = oth- 
erwise. The vF v spectral indices are denoted by a(> 0) 
and b(< 0). Note that S(x) can easily be generalized to 
accommodate spectral energy distributions (SEDs) with 
multiple components. 

The bolometric fluence of the model GRB for observers 
with 9 < 9j is given by 



F = 



dt 



At 



(3) 



where Xb is a bolometric correction to the peak measured 
vF v flux. If the SED is described by eq. ©, then X b = 
(a -1 — and is independent of e p k- The beaming- 

corrected 7-ray energy release £* 7 for a two-sided jet is 
given by 

p 

£* 7 = 47rd|(l - (ij) — 
where the luminosity distance 



+ z 



(4) 



, / s c . . T dz' 
d L ( z ) = —(l + z) / - 



(5) 



for a ACDM universe. Substituting eq. ([3]) for F into 
eq- © gives the peak flux 

£,.. 

fe P k = 



4nd 2 L (z)(l - Hj)At* Xb 



(6) 



Finally, by substituting eq. © into eq. ([TJ) , the energy 
flux becomes 



/.(*) 



£*-, Hjwfij,!) H{t;0,At)S(3 

A%d 2 rXz)(l 



(7) 



Because of the reduced spectral sensitivity of Swift to 
high-energy (> 200 keV) emission and the measurement 
of e p k in many GRBs, and to avoid complications from 
the potentially large number of parameters in our GRB 
model, we consider in this paper the simplest GRB SED 
with a = b = 0, so that Xb —>■ hi(x max /x m i n ) (see eqs. [3], 
[6] and [7]). The flat part of the spectrum takes the value 
of peak flux / £pfc given by eq. © , and is assumed to inter- 
sect the sensitive range of the GRB detector. By making 
the above approximation, we modify the fractional num- 
ber of bursts that occurs within a given redshift z and 
within a jet opening angle 6j (see eqs. [T5], [19] , and [2"T]); 
however, the physical interpretation of the problem is ex- 
pected not to change drastically when a, b ^ 0, which we 
plan to address in a subsequent paper. In the flat spec- 
trum approximation, a vF v spectrum (a = b = 0) has 
to be cut off at low (x m i n ) and high (x max ) energies in 
order to avoid divergent flux. The model spectrum is 
not likely to extend over more than two orders of mag- 
nitude, so that Afc < 5. Thus we take Ah, the bolometric 
correction factor, equal to 5 in our calculations. 

2.2. Bursting Rate of GRB Sources 

To calculate the redshift and the jet opening angle 
distributions, we need to determine the rate at which 



GRBs occur per steradian. For a ACDM cosmology with 
n m = 0.27, n A = 0.73, and H = 72 km s" 1 Mpc" 1 
(|Spergel et al.ll200"3t ). the event rate per unit time per sr 
for bursting sources is 



dN 

In 



c 

d£. 



^l( z ) ^co (£*7i Mji e pfe*j At*, z) 

, 7 d/ij d/i de p k* d(At*) dz 

(i + z) 3 x /o _ ;(i + z)3 + r! A ' 



(8) 



(derivation of this result is given in Appendix A). Here, 
n C o{£*~f , fJ>j > £pfc* , At* , z) — ?i co (z) Y(£*^ , fjj , c p k* , At* ) is 
the distribution of GRB sources differential in £* 7 , m, 
e p k*, At*, and z, and separability of the comoving rate 
density n co (z) — h co Y, SFR (z) of GRBs from their prop- 
erties is assumed, where S SFH (z) is GRB star forma- 
tion rate function. This is the crucial assumption of 
this treatment, and so needs restated. There is assumed 
to be no evolution of the emission properties of GRB 
with cosmic time. Thus, for example, the mean jet 
opening angles, the mean energies, and the mean val- 
ues of e v k* (whic h could be c orrelated with £*y to sa tisfy 
the lAmati et all (|2002f) and iGhirlanda et all (|2004fl re- 
lations when a, b ^ 0) do not change through cosmic 
time. The phenomenological validity of the the Amati 
and Ghirlanda relations for GRBs at different values of 
z suggests the plausibility of this assumption. These re- 
lations can be introduced in more detailed treatments 
by considering source spectral properties where e p k* is 
correlated with £* 7 . 

Employi ng an int egral formulation for beamed 
sources (|Dermerl Il992t 2006), the observed directional 
event rate for bursting sources 3 with vF v spectral flux 
greater than / ? at observed photon energy e is given by 



dN(> h) 



dn 


H Jf_ 








d/ij / d[i 




J in 




fc\Z> £*7 



de 



P k* l d(At*) 




dz 



d 2 L (z)n co (z)Y (£*^ , fij , e pk * , At*) 
(l + z) 3 y/ti m (l + z )* + nT~ 



(9) 

This expression can also be written in terms of a Heavi- 
side function involving a minimum flux threshold condi- 
tion f e > /g. For discrete values of £* 7 , e p k*, and At*, 

Y(£*j, fij , e p k*, At*) — g(fij) S(£*j — £* 1 q) 5{e p k* ~ e pfe*o) 
x 5{At* - At* ) , (10) 

where g(fJ-j) is the jet opening angle distribution, and 

fe is the instrument's detector sensitivity. The sensitiv- 
ity of a detector to a GRB depends on many factors, 
namely, the detector area, the detector efficiency, the 
fraction of the detector that is active, the fraction of 
the coded mask that is open, the average solid angle, 

3 In a more accurate treatment of detector response, one should 
calculate a photon-energy integration over the energy-dependent 
effective area, rather than describing a 7-ray telescope by a vF v 
flux sensitivity /g by the mean photon energy e of the sensitive 
waveband of the detector. A further complication, not considered 
here, is that fi may depend on the duration At of the GRB. 
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and th e internal background ([Band 2003) . iGehrels et al.1 
(2004) give the BAT's instrument detector sensitivity to 
be about ~ 10~ 8 ergs cm -2 s _1 , which is about an order 
of magnitude more sensitive than BATSE and Beppo- 
SAX, particularly in the case of those GRBs for which 
redshifts were obtained, as discussed below. 

We can make a rough estimate of the detector thresh- 
old for both Swift and pie-Swift instruments, noting that 
the mean flare size (cj>) — u(f> m - ul /(u — 1) for a power law 
size distribution with index u (> 1; u = 3/2 for a Eu- 
clidean distribution) . For Swift GRBs, we took the ratio 
of the measured GRB nuences to their (igo) durations 4 
to find that the mean flux of a detected Swift GRB is 
~ 5 x 10~ 8 ergs cm" 2 s -1 . For an index u = 1.1 - 1.2, 
as appropriate to the Swift range, we find that the Swift 
is sensitive to GRBs as weak as 10 -8 ergs cm" 2 s" 1 . We 
take this energy flux as the effective Swift flux thresh- 
old, recognizing that the actual detection of GRBs with 
Swift is mor e complicat ed, involving both a rate and im- 
age trigger (|Bandll2006f ). 

The assignment of an effective flux threshold for the 
pie- Swift detectors is complicated by the fact that they 
involve a variety of detectors — BATSE, Beppo-SAX, 
HETE-II, and INTEGRAL— each with its own specific 
triggering criteria. The mean flux of the pve-Swift sam- 
ple can be obta ined from the valu e s of fluence and 
t 90 compiled by iFriedman fe Blooml (|2005l ). We find 
that the mean flux of the pie-Swift GRBs is ~ 7 x 
10 -7 ergs cm -2 s _1 , meaning that the flux threshold sen- 
sitivity is a factor 5 less than th is value, because u = 
1.2 at fluxes > 10" 7 ergs cm" 2 s" 1 (Paci esas et alj [l999). 
For simplicity, we use the value of 10" 7 ergs cm" 2 s" 1 
for the mean sensitivities of the GRB detectors prior 
to Swift. Hence, in this work we assume the detector 
threshold for Swift and pie-Swift GRBs to be ~ 10" 8 
and ~ 10" 7 ergs cm" 2 s" 1 , respectively. 

Integrating over /i in eq. ((9]) gives the factor (1 — /ij) 
describing the rate reduction due to the finite jet opening 
angle. Hence, eq. ([9]) becomes 



dN(> h) 



dz 



d\(z) h co {z) 6(f' t 



fe) 



(1 + Zf y/SlmO- + Zf + r»A 



(11) 



where f e is given by eq. ©. 

In order to close the system of equations and solve 
for the redshift distribution, we must also specify the 
functional form of g(/ij) as well as the comoving GRBs 
rate density h co (z). GRBs with small opening angles, 
that is, with 8j <C 1, will radiate their available energy 
into a small cone, so that such GRBs are potentially 
detectable from larger distances (though not always, as 
a different portion of the spectrum is observed), with 
their rate reduced by the factor (1 — (ij). By contrast, 
GRB jets with large opening angles are more frequent but 
only detectable from comparatively small distances. The 
intrinsic opening angle distribution can be much different 
than the mea sured opening-angle distribution because of 
these effects (jGuetta. Piran. fe Waxmanll2005t ). 

4 swift.gsfc.nasa.gov 
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Fig. 3. — The number of sources per unit 9j angle for s = —1.2, 
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radians. 

Since the form for the jet opening angle g(fJ-j) is un- 
known, we consider the function 



a(P]) = .90 (1 - VjY H(nf,n }1 



mm? Mjmaxy j 



(12) 

where s is the jet opening angle power-law index; for a 
two-sided jet, //j m i n > 0. Because the distributions are 
normalized to unity, 



1 



9o 



(l /^jmin) (1 Mjmax) ~^~ & 



(13) 



Equation (| 1 2[) hardly exhausts the possible expressions 
for g([J,j), but contains interesting asymptotic behaviors. 
In Figure El we illustrate the nature of this function 
by fixing the range of the jet opening angle for differ- 
ent value of s. Here we take 9j lying between 0.05 and 
0.7 radians (w 2.8° - 40°). This figure shows that the 
fraction of bursts with small jet opening angles 9j <C 1 
increases with decreasing values of s < 0. Therefore 
most bursts have large jet opening angle when s > 0. 
Since beaming effects are only important for negative s, 
and because beaming effects are certainly important for 
GRBs, we only consider negative value of s in our anal- 
ysis. Furthermore, from equation (|13p it is clear that 
if s < — 1, it is necessary to have a minimum opening 
angle 9j = arccos /Xj > in order to avoid a divergent 
fraction of GRBs with small jet opening angles. Note 
also that g(/J>j), eq. fp~2j) . approaches a S- function when 

Mjmin ' Mjmax- 

It has bee n sug gested (e.g., iTotanil Il997t 
iNatarajan et al.l [2005) that the GRB formation 
history is expected to follow the cosmic SFR (S SFR ) 
derived from the blue and UV luminosity density of 
distant galaxies, though with differences related to the 
metallicity dependence and fraction of the subset of 
high-mass stars that are GRB progenitors. However, as 
we show in § 3, in order to obtain the best fit to the 
Swift redshift and the pie-Swift redshift and jet opening 
angle samples, we have to employ a GRB formation 
history that displays a monotonic increase in the SFR 
at high redshifts (SFR5 and SFR6; see Fig. gj). To 
investigate the nature of the GRB formation history, in 
this work we also assume that n co {z) oc E SFR (z), and 
consider specific functional forms for S SFR (2;). Due to 
the large range of uncertainties in the SFR at z > 1, 
we first consider three different functional forms for 
the SFR. For the lower SFR (LSFR) and upper SFR 
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(USFR) that should bound the actual SFR, we consider 
the analytic function 



l + ai 



(1 + z)- a * + ai(l + z) a 3 ' 



(14) 



with ai = 
a 3 = 3.0(3.0) for 



(jWick. Dermer fc Atovanl 120041 ) . 
0.005(0.0001), a 2 = 3.3(4.0), and 
the LSFR (USFR). The LSFR (SFR2) and USFR 
(SFR4) describe extreme ranges of optical/ UV measure- 
ment s without and with dust extinction (|Blain et al.l 
1999) corrections, respecti vely. We also cons i der th e 
SFR history (SFR3) of IHopkins fc Beacoml (2006), 
which is intermediate between the LSFR and USFR, 
using their analytic fitting profile given by 



t(*) = 



1 + (a 2 z/ai) 
1 + (z/a 3 ) a * 



(15) 



where a\ = 0.015, a 2 = 0.10, a 3 = 3.4, and = 5.5 
are their best fit parameters 5 . Figure 2] shows these 
SFR functions, with the function normalized to unity 
at the present epoch z = 0. The constant comoving 
density, SFR1, is also plotted. Note that SFR2, SFR3, 
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Fig. 4. — Star formation rate (SFR) function for GRBs, as- 
sumed to be proportional to comoving SFR histories as shown. 
The solid line (SFR1) is a constant comoving densi ty; SFR2 and 
SF R4 are the lower and upp er SFRs given by eq. (1 1 3 1) ; SF R3 is 
the IHopkins fc B cacom (2006) SFR history given by eq. 1141 1. with 
(al = 0.015, a2 = 0.1, a3 = 3.4, a4 = 5.5); and SFR5 (al = 0.015, 
a2 = 0.12, a3 = 3.0, a4 = 1.3) and SFR6 (al = 0.011, a2 = 0.12, 
a3 = 3.0, a4 = 0.5) are the SFRs that give a good fit to the 
Swift and pie- Swift redshift distribution. 



and SFR4 are similar when z < 1, but differ greatly at 
higher redshifts. We also consider two alternative func- 
tions shown in Fig. |H SFR 5 and SFR6, that display a 
monotonic increase in the SFR at high redshifts repre- 
senting positive evolution of the SFR history of the uni- 
verse to z > 5. Later in the analysis we will show that 
for self-consistency we have to utilize SFR5 and SFR6 
instead of SFR2, SFR3, and SFR4, to obtain the best fit 
to the Swift redshift, and the pre-Swift redshift and jet 
opening angle samples. 

2.3. Redshift, Size, and Jet Opening Angle Distributions 

5 Adjustments to these values in later versions of the paper by 
Hopkins & Bcacom (2006) produce negligible differences in our re- 
sults. 



The directional GRB rate per unit redshift with energy 
flux > fg is thus given by 



ego 



d 2 L (z) h co (z) 



dN(> h) = 

'/° dz H (2 + s) (1 + 2 )3 v /n m (i + z)3 + n A 

X {[l-max^Mjmin)] 2 ^- (l-ftmax) 2+S } , (16) 

after substituting eq. (I12p into eq. (1 1 1 1) and solving; here 

A* = 1 ^ ^- • (17) 

A-Kd\{z) At* f s Xb 

When f e — f g , where f g is the vF v flux threshold sensi- 
tivity of a GRB telescope that observes in a waveband 
centered at e — e, eq. (fl"6]) gives the redshift distribution. 
The size distribution of GRBs in terms of their vF v flux 
/ e is then simply 



dN(> A) 

dn 



ego 



H (2 



dz 



d\{z) h co (z) 



(1 + zf V n m(l + z) 3 + fiA 



X {[1 - max(A J ,/ijmin)] 2+S ~ (1 ~ Mjmax) 2+S } ■ (18) 

The maximum redshift z max can be obtained by solv- 
ing the luminosity distance di(z) using equation (|17p to 
obtain 

(19) 
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Finally, the jet opening angle distribution (see eq. [TT| ) 
is given by 

dN(> f t ) 



dfldfij 



Ho 



S'OjXl - Mi) 



(20) 



dz 



d\(z) n co (z) 



(1 + zf v/n m (l + z) 3 + OA ' 

where g(fJ-j) is given in equations (|12[) and (|13[) . and the 
upper limit of the redshift integration is determined by 
equation (fT9|) . 

3. RESULTS 

In our GRB model, there are seven adjustable pa- 
rameters: the vF v spectral power-law indices a and 6, 
the power-law index s of the jet opening-angle distribu- 
tion, the range of the jet opening angles #j m i n and f?j max , 
the absolute emitted 7-ray energy £* 7 , and the detec- 
tor threshold f g . As already mentioned, we consider, for 
simplicity, a flat vF v GRB SED, so that a = b = 0, leav- 
ing only the bolometric correction factor which is set 
equal to 5. The remaining parameters #j m i n , 6*j max , s, 
and £* 7 are constrained by observations, and the detec- 
tor properties are used to estimate / ? . 

Besides these parameters, we must also assume a GRB 
comoving rate density h co (z) — ti co Yisfr{z), which is 
also constrained by fitting our model to the observed 
data. Adopting the uniform jet model, we assume that 
the energy per solid angle is roughly constant within 
the well-defined je t opening angle (e.g.. lFrail et al.ll2001b 
iBloom et al.ll2003h . With these assumptions, we investi- 
gate the behavior of the redshift and opening angle dis- 
tributions. Before we proceed to compare our model of 
redshift and jet opening angle distributions with the sam- 
ples obtained by the Swift and prc-Swift instruments, we 
need to explore some of the above parameter space. 
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Fig. 5. — Directional event rate per unit redshift [one event 
per (10 28 cm) 3 per day per sr per z] for an assumed 5-function 
distribution of the jet opening angles, so that g(fij) = — 



jij), with 9j 



= 0.1 radian (solid-curve), using SFR3 and /j 
-2 „-l 



with #j max = 0.1 radian, and 0j m i n 
with s = 0. 



0.05, 0.07, and 0.09 radians, 



3.1. Parameter Study of GRB Model 

For illustrative purp oses, we utilize SFR3 
([Hopkins fc Beacoml [2006) for the GRB comoving 
rate density. Furthermore, we also assume that s = 0, 
#jmin = fjmax = Oj =0.1 radian (this assumption implies 
g(fij) is a (5-function, see eqs. [II] . [12] , and Fig. [5]), 



At, = 10 s, £„ 7 = 2 x 10 51 ergs, and /<? = 10" 8 and 
10~ 7 ergs cm~ 2 s -1 for the flux sensitivity of detectors 
for Swift and before Swift, respectively. We choose 
At, = 10 s, and £* 7 = 2 x 10 51 ergs because these are 
the mean values of the GRB duration measured with 
BATSE assuming that BATSE GRBs are typically at 
z w 1, and the mean beaming corrected 7-ray energy 
release (see Fig. 1), respectively. We shall use At* = 10 s 
throughout this work. The calculation is normalized 
to a local [z <C 1) event rate of one GRB per day per 
(10 28 cm) 3 , so that the normalization of the local rate 
density is h co = 1.157 x 10~ 89 cm~ 3 s _1 . 

With the jet opening angle 8j — 0.05 and 8j =0.1 radi- 
ans, and the detector threshold fg = 10 -7 ergs cm~ 2 s _1 , 
the maximum redshifts from which model GRBs could 
be detected are z max = 5.4 and 3.02, respectivley (see 
eq. \19\). The directional event rate per unit redshift 
is shown in Figure Setting #j max = 0.1 radian, and 
#jmin = 0.05, 0.07, 0.09 radians (the three dashed curves 
in Figure [5]) , we plot the directional event rate per unit 
redshift for different values of 0j m [ n (see eq. [TH]). As 
$jmin - > 0jmax, the directional event rate per unit redshift 
of the respective #j m i n -curves approaches the directional 
event rate per redshift of the 5-function g(fJ.j) with 9j = 
0.1 radian. In Figure [5] we examine three different val- 
ues of the jet opening angle, 9j mm — #j max = 9j = 0.05, 
0.1, and 0.2 radians, with two values of flux sensitivity, 
namely 10 -8 and 10~ 7 ergs cm~ 2 s _1 . As can be seen, 
the smaller the jet opening angle the larger the maxi- 
mum redshift from which the GRB is detectable. The 
maximum observable redshift z associated with a par- 
ticular opening angle is indicated by the vertical line in 
each curve. Moreover, Figure [5] also shows that what- 
ever the lower sensitivity telescope detects will also be 
detected by a telescope with better sensitivity, because a 
flat vF v spectrum is considered. Thus Swift, being more 



140 

3 20 

„ 100 
o 

x 80 

CI 
TJ 

n 60 

T3 
•g 

■O 40 
20 




0; = 0.2 rad 



0.1 



/ 



0.05 



1 2 3 4 5 6 

Redshift z 

Fig. 6. — Directional event rate per unit redshift, in units of 
GRBs per s per sr per unit redshift z, normalized to one event 
per (10 2s cm) 3 per day, with 9j = 0.05, 0.1, and 0.2 radian, 
using SFR3. The thin and thick curves represent the instru- 
ment's energy flux sensitivity of 10 -8 ergs cm~ 2 s _1 (Swift) and 
10 — 7 ergs cm -2 s — 1 (pre-Swift), respectively. 

sensitive than earlier telescopes, would in this model de- 
tect a sample of GRBs encompassing those detected with 
instruments before Swift. Spectral behavior can be con- 
sidered in more detailed treatments and will change the 
relations for maximum distances, though generally with 
better sensitivity the telescope can view to larger dis- 
tances, especially when the more sensitive telescope can 
observe lower energy photons. 



s = 0.0 




Fig. 7. — Directional event rate per unit redshift [one event 
per (10 28 cm) 3 per day per sr per z] with 9\ m i n = 0.05 rad and 

-0.9, -1.5, 



: 0.7 rad, 



0.0, -0.5, 



jmm 

2.5, and -10.0, with 



SFR3 and /j = 10 -8 ergs cm" 2 s" 



By varying the opening angle power-law index pa- 
rameter s over values s = 0.0, —0.5,-0.9, —1.5, —2.5, 
and —10.0, with 6*j m j n = 0.05, 6j max — 0.7 radians and 
ft = 10~ 8 ergs cm~ 2 s _1 , while keeping all other assumed 
parameters as before, Figure [7] shows that the burst rate 
per unit redshift is progressively distributed to higher z, 
as expected, since the more negative s is, the larger is 
the fractional number of bursts that occurs at small jet- 
opening angles 8j (see also Figured]). More importantly, 
as s <gC 0, the directional event rate per unit redshift 
reaches an asymptotic limit (see plots for s = —2.5 and 
-10.0 in Figure [7]). 

The final parameter that we need to explore is the 
corrected 7-ray energy £* 7 . By setting £* 7 = 5 x 10 50 , 
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Fig. 8. — Directional event rate per unit redshift [one event 
per (10 28 cm) 3 per day per sr per z] with 0j m i n = 0.05 rad and 
6» jmax = 0.7 rad for £„ 7 = 5 X 10 50 , 3 X 10 51 , and 10 52 ergs, using 
SFR3. 

2 x 10 51 , 4 x 10 51 , and 10 52 ergs with 9 jmin = 0.05 and 

0jmax = 0.7 radians, s = 0, and f s — 10~ 8 ergs cm -2 s -1 , 
Figure [5] shows that as the emitted 7-ray energy in- 
creases, the bursting rate is shifted to higher z, as ex- 
pected. 

3.2. Parameter Fits to the Swift and pre-Swift 
Redshift samples 

From iFriedman fc Blooml (|2005h . the pre- Swift jet 
opening angle distribution extends from 6j = 0.05 to 
9j = 0.6 radians. However, other workers have reported 
that the observed je t-opening angles are as large a s 
6j = 0.7 radians (e.g.. iGuetta. Piran. fc Waxmanl 12005). 
Hence, in this work we explore different ranges of 9j, 
ranging from 9j m - m = 0.05 radians to 9- jmax = tt/2, 0.7, 
and 0.4 radians. We also use At* = 10 s as the average 
GRB durati on in the stationary fram e. We first con- 
sider SFR3 ([Hopkins fc BeacomlbooH to fit the redshift 
distribution from the Swift and pie-Swift samples, and 
the opening angle distribution from the pie-Swift sam- 
ple. A large number of jet opening angles associated with 
the Swift redshift sample is lacking, so our calculated jet 
opening angle distribution will be used as a prediction for 
Swift. As indicated in Figure [TJ the corrected 7-ray en- 
ergy release is broadly distributed around the mean value 
of £* 7 = 2 x 10 51 ergs; hence we can explore different 
values of 7-ray energy that best fit the observed redshift 
and opening angle samples. It is important, however, 
that the best-fit value is not too far off from the mean 
value in order for our model to be in accord with obser- 
vations. Later in the analysis we will show that our best 
fit corrected 7-ray energy £* 7 is indeed only a factor of 
2 larger than the mean. 

In Figures [9[ [10] and [Til we plot the calculated red- 
shift distribution at three different corrected total 7-ray 
energies, £* 7 = 2 x 10 51 , 5 x 10 51 , and 10 52 ergs, respec- 
tively. In each panel, each curve represents the calcu- 
lated cumulative redshift distribution for the power-law 
index s of the jet opening angle distribution. The values 
s = —0.6, —1.2, —1.8, and —2.4 are plotted, from far 
left to far right, respectively. In each figure, from top to 
bottom, we have #j m in = 0.05 and 9- irn3X = ir/2 radians, 
#jmin = 0.05 and 9j max — 0.7 radians, and 9 ]min = 0.05 
and #j max = 0.4 radians, respectively, as the jet opening 



angles. The cumulative redshift distribution (bold solid 
curves) of the Swift and pie-Swift samples are plotted 
for comparison. We do not plot s < —2.4, since there 
is little variation in the redshift distribution (see Fig. [7]). 
The results in Figures (H [UJl and QT] clearly show that 
the redshift distributions for different values of s from 
the model do not fit the Swift redshift sample. Further- 
more, the fits to both the Swift and pre- Swift samples do 
not correlate together with the associated range of the jet 
opening angle. In fact it indicates that we need larger 
burst rates at higher redshifts (z > 3) . These results rule 
out the lower SFR, SFR2, as can be seen in Figure [P2l 

It is important to note at this point that the fit to the 
pie-Swift and Swift redshift distribution should rely on 
a single comoving rate density of GRB with the same 
physical parameters. If one particular GRB rate model 
cannot fit the observed pre-Swift and Swift redshift dis- 
tribution at the same time, then that GRB model has to 
be eliminated, unless additional complications are intro- 
duced into the model, such as the change of GRB model 
parameters with time. 

The comparisons also rule out the possibility that the 
upper SFR explains the GRB data, even though SFR4 
has a higher SFR than SFR2 at large redshifts. In Fig- 
ure[l3]we fit the pre-Swift and Swift redshift samples us- 
ing three different 7-ray energies, namely £* 7 = 5 x 10 50 , 
2 x 10 51 , and 10 52 ergs, with s = -0.6, -1.2, -1.8, and 
—2.4 using SFR4. The results for the range of the jet 
opening angle #j m i n = 0.05, and 9j max = 0.7 radians 
show that the calculated redshift distribution between 
Swift and pre-Swift are anticorrelated, that is, smaller 
values of s, which would improve the fit for Swift GRBs, 
makes the fit for pre-Swift GRBs worse. However, with 
s = —0.6 and £„ 7 = 10 52 ergs (two bottom panels), the 
redshift distributions are marginally acceptable by the 
model, noting that the Kolmogorov-Smirnov one-sided 
test statistic for a sample of 41 and 31 GRBs at the 10% 
level are 0.167 and 0.187, respectively (that is, only 10% 
of randomly chosen samples should exhibit maximum 
vertical deviations from the model cumulative distribu- 
tion; if it is greater than these values, then the model 
is an inadequate representation of the data at the 90% 
confidence level). Only 2.5% of randomly chosen sam- 
ples should have vertical excursions greater than 0.212 
and 0.238 for a sample of 41 and 31 GRBs, respectively. 

When we calculate the opening angle distribution to 
fit the pre-Swift opening angle sample in Figure I14| 
however, the fit is adequate for £» 7 = 10 52 ergs and 
s = —1.2 but not for s — —0.6. Hence, the upper SFR 
(SFR4) is statistically eliminated. Nevertheless, the re- 
sults in Figure [T3l suggest that a greater GRB burst rate 
at larger redshifts is required to fit the redshift distri- 
butions. The cumulative jet opening angle distribution 
(bold solid curves) of the pre-Swift sample is plotted for 
comparison, but not for Swift, since the available open- 
ing angles for the associated Swift redshift sample are 
lacking. 

As a result, we consider a GRB comoving source rate 
density that exhibits positive evolution above redshift 
z > 3. We therefore modify the Hopkins & Beacom 
SFR (SFR3) to obtain the comoving rate-density that 
attains (1) a constant comoving rate-density at high red- 
shift z (SFR5), and (2) a monotonically increasing co- 
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moving rate-density with increasing redshift (SFR6) in 
order to provide a high GRB burst rate at large z (see 
caption to Figured] for coefficients of the SFR models). 

In Figures[T5]and[T6]we begin to see some improvement 
in fitting the pie- Swift and Swift redshift samples when 
using SFR5 and SFR6 with £„ 7 = 2 x 1 51 ergs. Never- 
theless, the fitting is still somewhat inconsistent between 
the pie-Swift redshift distribution and the jet opening 
angle. For example, in Figure [TBI using SFR5, the fitting 
to the pre-Swift and Swift redshift distribution shows a 
good fit with s = —1.2 for a range of jet opening an- 



gles between ft^ 



0.05 and ft. 



0.4 radians (two 



bottom panels). However, the result from the fitting be- 
tween the pre- Swift redshift distribution and the opening 
angle distribution with s = —1.2 is inconsistent (see Fig- 
ures [TS] and rTB]) . Above redshift z > 2, the calculated 
redshift distribution underestimates the fractional num- 
ber of bursts that should occur at high redshifts z. On 
the other hand, the calculated opening angle distribution 
overestimates the fractional number of GRBs that should 
occur at small opening angles. Recall that if a burst oc- 
curs at a small opening angle, then this translates into 
the possibility of observing this particular burst at a high 
redshift z. Hence, the fit with the range of the jet open- 
ing angle 9j m in = 0.05, and 0j max = 0.4 radians is not yet 
acceptable. 

In Figures [T7] and [18] we see similar results when us- 
ing SFR6. However, it is interesting to note that we 
can resolve this problem if we could increase the frac- 
tional number of GRBs at both high and low redshifts 
z. As a result, this will provide a self-consistent fitting 
between the pie-Swift redshift and the jet opening angle 
distributions, as well as the Swift redshift distribution. 
This problem can be solved by adjusting the 7-ray en- 
ergy £* 7 . Using £* 7 = 4 x 10 51 ergs, Figures [151 and 
120] and Figures [2T1 and 1221 show a statistically acceptable 
fit to the pre-Swift and Swift redshift distribution and 
the pre-Swift opening angle distribution with s = —1.2, 
and the range of the jet opening angle #j m in = 0.05 and 
ftjmax = 0.7 radians for SFR5 and SFR6, respectively. 
For SFR5, the fitting can be improved with s = —1.3, as 
indicated by the solid curves (see Figures [19] and [20]) . 

The above analyses show that £* 7 = 4 x 10 51 ergs gives 
the best fit to the pre-Swift redshift and the jet open- 
ing angle samples, and the Swift redshift sa mple. This 
best-fi tted 7-ray energy is int eresting since iFrail et all 
(|2001| ) and iBloom et all (|2003f ) show that the standard- 
energy reservoir is £* 7 = 10 51 er gs. However, with 



a larg er sample set, the data from iFriedman & Bloom 
(2005) shows that £» 7 is broadly distributed between 



£* 7 ~ 1Q51-52 er g Sj as shown in Figure [T] Because we 
have been able to obtain an acceptable fit to the data 
with discrete values of £ * 7 and At* , a more complicated 
model where a range of values of £* 7 and Ai» were al- 
lowed would certainly permit excellent fits to the data. 



3.3. Estimating the Beaming Factor and 
the Luminosity Function 

We can utilize our best-fitted parameters to estimate 
the average beaming factor < f^ 1 >. The beaming fac- 
tor average over the opening angle distribution is given 



by 



<f b > 



(21) 



2 + s 



2+sl 



(1 Mjmin) (1 /^jmax) 

1 - //■ ■ ) 1+s - (\ - U ' ) 1+s 



where g(fij) is given in equation (fl2]) . Using the best- 
fitted parameters s = —1.2 and s = —1.3 and the range of 



the jet opening angles between #j, 



0.05 and ft. 



0.7, we obtain < /& >~ 1 ~ 34 — 42 as the beaming factor 
for s — —1.2 and s = —1.3, respectively. Our value of 
< fb > _1 is about a factor of 2 smaller than the value 



obtained bv lGuetta. Piran. k, Waxmanl (|2005f ). 

We can also infer the luminosity function from the jet 
opening angle distribution. From equations {H} and ([9]), 
the jet opening angle distribution is dN/dfij = (?(/%), 
which is related to luminosity by the expression 



dN 
dLZ 



dfij 



dL* 



(22) 



where L» is the GRB apparent isotropic luminosity, 
and g(fJ-j) is the jet opening angle distribution given 
in equation (|12[) . Since the energy flux is defined as 
$_e = L*/47T<i| or f t = vF v = v*L VJf / '47rd|, then the 
isotropic luminosity in the bursting frame is = is*L u *, 
so that = 4nd 2 L f e . Utilizing equation (TIT]) , the jet 
opening angle in terms of isotropic luminosity is given 

by 

f 



Mi = 1 - 



(23) 



Substituting equation ([23]) into equation ([22]) . we ob- 
tain dN/dL* cx ( s+2 ) after utilizing the jet open- 
ing angle distribution in equation (|12[) . Since our 
best fit to the pre-Swift redshift and opening angle 
samples and the Swift redshift sample constrains the 
jet opening angle distribution power-law index s be- 
tween — 1.3 and —1.2 for using SFR5 and SFR6, re- 
spectively, our apparent isotropic luminosity function is 
cx L~ 3 ' 25 within the uniform jet model. In the uni- 
versal structured jet model, the luminosity function is 
sugge s ted to be cx L~ 2 (e.g . . iGuetta. Piran. fc Waxmanl 
I2005t iPerna. Sari. & Frail l2003f) : others have sug- 
gested that the luminosity function is a L ~ 23 (e.g., 
iLloyd-Ronning. Frver. fc Ramirez-Ruizl [2002) . For our 
best fitting model with jet opening angles ranging from 
8j = 0.05 to 8j — 0.7 rad, this implies bolometric 7-ray 
luminosities in the range w 2 x 10 51 to w 3 x 10 53 ergs 
s" 1 . 

The results for our best fitting model are shown in 
Figure [23] The fraction of high-redshift GRBs in the 
high-redshift universe implied by our fits using SFR5 and 
SFR6 is estimated at 8 - 12% and 2.5 - 6% at z > 5 and 
z > 7, respectively. If this fraction is not observed with 
Swift, then it suggests that SFR5 and SFR6 do not con- 
tinue with positive evolution z>5. Indeed, it is likely 
that the formation rate of GRBs at high redshift declines 
at sufficiently high redshift. We can also derive the frac- 
tional number of GRBs with z < 0.25 from Figure [23] 
We find that < 0.2% of long duration GRBs should occur 
in the low-redshift (z < 0.25) universe. 
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4. DISCUSSION AND CONCLUSIONS 

In this work, we considered whether the differences be- 
tween the pre- Swift and Swift redshift distributions can 
be explained with a physical model for GRBs that takes 
into account the different flux thresholds of GRB detec- 
tors. The model presented here parameterizes the jet 
opening angle distribution for an assumed flat vF v spec- 
trum, and finds best fit values for the 7-ray energy release 
for different functional forms of the comoving rate den- 
sity of GRBs, assuming that the properties of GRBs do 
not change with time. Adopting the uniform jet model, 
we assumed that the energy per solid angle is roughly 
constant within a well-defined jet opening angle. The 
pre- Swift redshift sample suggests that the range of the 
jet opening angle is between 0.05 and 0.7 radians. We 
explored three different ranges of jet opening angle with 
#jmin = 0.05, and 0j max = tt/2, 0.7, and 0.4 radians. 

The results in § 3 show that an intrinsic distribution in 
the jet opening angles and single values for all other pa- 
rameters yields a good fit to the pre-Swift and Swift red- 
shift samples, and furthermore provides an acceptable 
fit to the distribution of opening angles measured with 
pie-Swift GRB detectors. A good fit was only possible, 
however, by modifying the Hopkins & Beacom SFR to 
provide positive evolution of the SFR history of GRBs 
to high redshifts (SFR5 and SFR6; see Fig. gj). The 
best fit values were obtained with £* 7 = 4x I0 51 ergs, 
#jmin = 0.05, 0jmax = 0.7 radians, and s = —1.25. The 
best-fit value for the absolute 7-ray energy released by 
a GRB is a bout a factor of 2 l arger than the mean 
value of the iFriedman fe Blooml (|2005f ) sample of pre- 
Swift GRBs. 

The results of our fitting therefore indicate that GRB 
activity was greater in the past and is not simply pro- 
portional to the bulk of the star formation as traced 
by the blue and UV luminosity dens ity of the universe. 
This i s contrary to the conclusion of iPorciani fc Madaul 
(12001 based upon a study of BATSE size distribution 
and a small number of GRBs with known redshifts in 
the pre-Swift era. They obtain worse fits for SFRs 
that increase a t high redshift. The BATSE analysis of 
iMeszaros et al.l ()2006t ) gives results that are consistent 
with the req uired monoton i cally increasing GRB SFR 
that we find. iDaigne et alj (|2006f ) also conclude that a 
rising GRB SFR is implied, and also note that the star 
formation rate cannot increase in proportion to greater 
past GRB activity, as this would result in overproduc- 
tion of metals. Changes in the progenitor population or 
GRB properties, for example, due to changing metallic- 
ities with time, might account for the large number of 
high-redshift GRBs observed with Swift, but this would 
introduce a large number of unconstrained parameters. 
Moreover, one would then have to explain why the Amati 
and Ghirlanda correlations relating spectral properties 
and energy releases are insensitive to redshift. 

Assuming that the GRB properties do not change with 
time, we therefore find that the GRB source rate density 
must display positive evolution to at least z > 5. Why 
would GRB activity be greater in t he past? One pos- 
sibility suggested by recent studies (jStanek et al.l [20061 : 
iFruchteret al. 2006) on the metallicities of host galaxies 
of GRBs is that GRBs derive from progenitor stars with 
a lower metallicity than the progenitors of Type II SNe. 



In this case, one might expect that the peak of GRB 
activity would occur earlier in the universe than during 
the peak of star formation activity at z f» 2 - 3. Never- 
theless, the high-redshift fraction of GRBs detected with 
Swift will test the validity of SFR5 and SFR6 used to 
model the GRB data. 

The models that fit the data imply that 8 - 12% of 
Swift GRBs occur at z > 5, in accord with the data as 
shown in Figure |2"51 Our model predicts that 2.5 - 6% 
of GRBs should be detected from z > 7. If no GRB 
with z > 7 is detected by the time « 100 Swift GRBs 
with measured redshifts are found, then we may conclude 
that the GRB formation rate has begun to decline above 
z rs 5. By contrast, if a few z > 7 GRBs have been de- 
tected within the next year or two with Swift, then this 
would be in acco rd with our mod e l and would support 
the conjecture of iBromm fc Loe b (2006) that a second 
episode of Population I and II star formation took place 
at z > 5, possibly even including some GRBs that orig- 
inate from Population III stars. Our model parameters 
suggest that GRBs can be detected with Swift to a max- 
imum redshift z rs 20, even given the limits placed on 
the range of jet angles, in particular, 6 m i n > 0.05. 

The relationship between the galaxies that host GRBs 
and their metallicities remains controversial, with evi- 
denc e for high-redshift, low -metallicity GRB host galax- 
ies (|Fruchter et al.l l2006h countered by examples of 
GRBs found in galaxies with moder at e metallicity 
(jBerger et all 1200a iFynbo et al.ll2006ah . IStanek et all 
(2006) find that 5 GRBs, or several per cent of the 
GRB population with known redshifts, are found at low 
(z < 0.25) redshift in galaxies that have very low metal- 
licity compared with Solar metallicity and with the dis- 
tribution of metallicities of low-redshift galaxies in the 
Sloan Digital Sky Survey. In 4 out of the 5 low-z GRBs, 
the apparent isotropic energy release of these GRBs is 
much smaller than 10 51 ergs, and so these GRBs are 
unrepresentative of typical long duration GRBs. More- 
over, the large number of these low redshift GRBs al- 
ready means that they are not typical of the sample of 
GRBs considered in this paper, as the model prediction 
is that « 0.14 - 0.2% of the "classical" GRBs should 
be detected at low redshifts. 6 These low-redshift GRBs 
may also belong to a completely distinct population of 
GRBs compared to their high-redshift counterparts, as 
also suggested by th e differing lumino sity functions of 
the two populations (jLiang et al.ll2006l ). Alternatively, 
the low-redshift GRBs could be normal GRBs that ap- 
pear weak because, fo r example, their jetted emission is 
observed off-axis (e.g. lYamazaki et~aT1l2003[ ). 

Although we have restricted our treatment to the 
long-duration, soft-spectrum class of GRBs, the sepa- 
ration between the short hard and long soft class of 
GRBs is not distinct. Moreover, recent observations 
of nearby GRBs show that the nearby (z = 0.125) 
GRB 60614, a long durati on GRB in terms of its light 
curve ( Gehrel s et al.l I2006D . lacks a supe rnova excess 
(|Fvnbo et all 120061 : iDella Valle et al J 12009) . does not fit 
neatly into e ither the short hard or long soft classes 
(|Zhangl r2006h . and could represent a separate popula- 

6 A value of ss 0.3% can be estimated by considering the frac- 
tional volume within 2 = 0.25 compared with 2 = 1, and that the 
SFR activity at z = 1 was a factor = 5 greater than at 2 < 0.25. 
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tion of GRBs (iGal-Yam et al.l [20061. Phenomenological 
studies have also indicated that the long GRB popula- 
tion is bimodally distributed in terms of the autocorrela - 
tion fu n ction of the light curves (e.g., lBorgono"voll2004[ ) . 
iTavanil (jl998h argues for multiple populations by exam- 
ining the size distributions of different groups of BATSE 
GRBS separated according to spe ctral properties . From 
an analysis of durations of BATSE. [Horvath et alj (|2006l ) 
identifies an intermediate duration population that dis- 
plays the softest spectra and comprises « 11% of the 
total BATSE GRBs. 

That we were able to obtain good fits to the different 
data sets using a single GRB model may indicate that the 
contribution of sub-populations to the total number of 
GRBs is small. Alternately, the model parameters might 
have been adjusted to take into account the aggregate 
long duration GRBs. This can be tested in future work 
using a larger Swift data set that allows subgroups to be 
separately studied based on spectra or durations. 

Figure l2~4"a shows the model integral size distribution 
(see eq. [TB]) of GRBs predicted by our best fit model. 
The plots are normalized to the current total number 
of observed GRBs per year from BATSE, which is 550 
bursts per 4-zr sr exceeding a peak flux of 0.3 photons 
cm" 2 s" 1 in the AE = 50 - 300 keV band for the 
At = 1.024 s trigger time (jBandl l2002t ). From the 
model size distribution, we find that = 340 to 360 GRBs 
per year should be detected with a BATSE-type detec- 
tor over the full sky above an energy flux threshold of 
~ 10" 7 ergs cm" 2 s" 1 , or photon number threshold of 
0.625 ph cm" 2 s" 1 . This range of values is determined 
by the two model fits with SFR6 and SFR5, respectively. 
We estimate from our fits that ~ 1190 - 1370 GRBs take 
place per year per Air sr with a flux > 10 -8 ergs cm" 2 s , 
or 0.0625 ph cm" 2 s _1 . The field of view of the BAT 
instrument on Swift is 1.4 sr (jGehrels et al.ll2()04h . im- 
plying that Swift should detect w 130 - 150 GRBs per 
year; currently Swift observes about 100 GRBs per year. 
This minor discrepancy may be a consequence of using a 
flat vF u GRB spectrum in our model. 

We also show the differential size distribution 
of BATSE GRBs f rom the Fourth BATSE catalog 
(|Paciesas et al.|[T999() in comparison with our model pre- 
diction in Figure l2~4"b. As can be seen, our model gives a 
good representation of the size distribution of the BATSE 
GRB distribution within the statistical error bars, except 
for a slight over prediction of the number of the brightest 
GRB s (see also iDaigne et al.l 120061 : iBottcher fe Dermerl 
2000). These brightest GRBs typically originate from 
z < 1, possibly suggesting a slight reduction in the num- 
ber of GRBs with z <C 1 compared to the GRB rate 
densities considered. Below a photon number threshold 
of 0.3 ph cm" 2 s" 1 in the 50 - 300 keV band, the ob- 
served number of GRBs falls rapidly due to the sharp 
decline in t he BATSE trigger ef ficiency at these photon 
fluxes (see iPaciesas et al.l ll999). The size distribution 
of the Swift GRBs will extend to much lower values, 
~ 0.0625 ph cm" 2 s _1 , and we can use our model to pre- 
dict the peak flux size distribution to fit the Swift data, 
noting however that the Swift triggering criteria are more 
complicate d than a sim ple rate trigger, particularly near 
threshold (|Bandll2006| ). 

We can derive the burst rate p in our local universe 
by normalizing our best fit models to the BATSE re- 



sults of 340 - 360 GRBs per year full sky with peak 
fluxes exceeding 10 -7 ergs cm" 2 s" 1 . Comparing this 
value with the rate calculated from eq. (fT8"]) for our 
best fit models, we obtain local GRB rate densities of 
?^GRB = 9.6 Gpc -3 yr" 1 for SFR5 and t^grb = 7.5 
Gpc -3 yr" 1 for SFR6. Since the mean volume occupied 
by an L* galaxy in t he local universe i s w 200 Mpc -3 
([Lovedav et aljfl992t IWiiers et ailll998h . we obtain the 
local event rate in our Galaxy, assumed to be a typ- 
ical L* galaxy, to be p « 1.9 x 10" 6 and 1.5 x 10" 6 
events per year per L„ galaxy, for SFR5 and SFR6, re- 
spectively. This implies that about 1 event occurs ev- 
ery 0.6 Myrs in our Galaxy. Comparing this to the rate 
of Type Ib/c supernovae of ~ 1 every 360 years in our 
Galaxy (Cappellaro et al. 1999), we estimate that the 
rate of GRBs is only « 0.06% of the r ate o f SN Ib/c 
supernovae in the Galaxy. iBerger et alj {2003) find that 
< 3% of GRBs could originate from SN Ib/c by com- 
paring the radio brightnesses of SN Ib/c with that of 
SN 1998bw associated with GRB 980425. Our results 
fall comfortably within the observational limits. We also 
find that about 1 out of every 20,000 SNe make a GRB, 
co mpared to the value of 1 out of every 10 5 - 10 6 found 
bv IDaigne et~aTl (|2006D . 

If GRBs do occur at the rate of more than once per Myr 
in the Milky Way, depend ing on uncertain metallicity ef- 
fects (jStanek et al.l 120061 ) . then the Earth will intercept 
the beam of a GRB about once every 25 Myrs, given the 
beaming factor (f^ 1 ) ~ 40 obtained in our analysis. If 
a Milky Way GRB is located at an average distance of 
10 kpc, then the average fluence deposited by the GRB 
is on the order of w 10 7 ergs cm -2 . Strong astrobio- 
logical effects begin to be important when the GRB flu- 
ence is « 10 8 - 10 9 ergs cm" 2 , or when the GRB occurs 
withi n ss 1 - 2 kpc with its b eam pointing in our direc- 
tion (|Dermer fc Holmeal2005l ). This will happen perhaps 
a few times per Gyr, so that GRBs could have signifi- 
cant effects on terr estrial evolution (M elott et alj 120041 
iThomas et aTll2005l ). 

The number of jet opening angles measured with 
follow-up observations of Swift GRBs is not yet suffi- 
cient to produce a statistically reliable sample. Thus our 
jet opening angle distribution shown in Figures [20] and 
[22] with s = -1.3 and -1.2, for SFR5 and SFR6, respec- 
tively, provide a prediction for Swift. With Swift and pre- 
Swift thresholds of ss 10" 8 and « 10" 7 ergs cm -2 s _1 , 
our model predicts the mean measured jet opening an- 
gle {Of) w 0.16, and 0.12 rad, respectively, which 
can be compared with the measured average value of 
(6j) « 0.125 for the pie-Swift sample. Thus we expect 
to detect more faint low-redshift, large opening angle 
GRBs that pie-Swift sate llites could not detect (see also 
iPerna. Sari, fc Frailll2003l ). It is interesting to note that 
with a better detector sensitivity, the mean jet angle is 
shifted to a higher angle and not vice versa, as might 
be expected. This suggests that on average we expect 
to see more long-duration bursts with larger opening an- 
gle (9j > 0.12 radian), meaning that the mean times for 
the achromatic breaks in the light curves could be longer 
than for the pie- Swift GRBs when one also takes into 
account the larger mean redshift of the Swift sample. 

In conclusion, we have developed a physical model for 
GRBs and obtained parameters to the model by fitting 
the GRB redshift and opening angle distributions. We 
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could only obtain statistically acceptable joint fits if the 
comoving rate density of GRBs increases monotonically 
to z > 5. Our fitting results give a rate of w 1 GRB 
every 600,000 yrs in our Galaxy if metallicity effects do 
not play a large role in the galaxies in which GRBs are 
formed. In this case, the increase in the GRB rate den- 
sity to high redshifts would have a different cause, for 
example, an episode of enhanced star formation at z > 5 
(jBromm fc Loebll2006D . The opening angle distribution 
of GRBs measured with Swift, and the fraction of low-z 



and high-z standard long-duration GRBs will test this 
model. 
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APPENDIX 

BURSTING RATE OF BEAMED COSMIC SOURCES IN A FLAT ACDM UNIVERSE 

The differential number of bursts per unit proper time per unit volume with beaming-corrected absolute 7-ray energy 
release between £, 7 and £» 7 + d£* 7 , jet opening cosine angle fij = cos8j between fij and fj,j+dfij, peak energy between 
e p k* and e p k* + de p k*, burst duration between At* and At, + d(Ai*), that occur in volume element between V* and 
V* + dV* within observer angle \jl = cos 9 between [i and [i + d\i during time interval between t, and t* + dt* is given 

by 

dN = h*(a,t*)dadndV*dt* . (Al) 

Here the parameter a represents the various parameters of the problem, namely £* 7 , fj,j, e p k*, and At*. Hence, the 
directional burst rate per unit observer time t is given by 



dN n*(at*)dad^dV* 

~dri ~ dn 



dt* 



dt 



(A2) 



In a flat ACDM universe, the Robertson- Walker metric describes the line element ds for a homogeneous, isotropic 
expanding universe, and is given by 



ds z =c z dt l -R\t) 



1 — kx 2 



(A3) 



where R(t) is the dimensionlcss scale factor of the universe and (x, 9, 4>) are the comoving coordinates of a point at 
rest in the Hubble flow. The curvature index k — —1/1Z 2 , 0, and +1/1Z 2 represents an open, flat, and closed universe 
model, respectively, where 1Z is the radius of curvature at the present epoch. At the present epoch, we let R(t) = 1 
and x correspond to the physical distance between two comoving points. The ratio of the scale factor, i?(t*), at the 
emitted time, t*, to the scale factor R — R(t ), at the present time, to, is given by, 

i?* 1 A, 



R 1 + z 

The following relations hold from the above definition, namely 

v* At 
v At* 



1 + z 



(A4) 



where the last expression applies to photons and relativistic particles in an expanding universe with R* = i?(t*), and 
R = R{t ). 

To calculate the differential volume element dV, in equation (|A2j) , we utilize the fact that the product of differential 
length elements corresponding to infinitesimal changes in the coordinates x, 9, and <f> equals the differential volume 
element dV . Hence, the proper volume element of a slice of the universe at time f* is, from eq. (|A3|) . 



dV* 



(R*x) 2 R*dxdQ* 
VI - kx 2 



Rlx 2 dxd£l* , 



(A5) 



specializing to a flat universe, where df2* represents the solid angle at the emitted time. Substituting equation (y\ > 
into equation (|A2|) for dV*, we obtain 



— — = n*{a,t*)R^x 2 dxdad^x 



dt* 



dt 



(A6) 



where in a homogeneous isotropic expanding universe with no intervening gravitating matter to distort the paths of 
the light rays, dfl* = dfl. 

Light travels along a null geodesic, ds = 0, which is a path along which 9 and 4> are constants. Thus in the Robertson- 
Walker metric the proper time element elapsed for a photon to propagate a differential comoving distance at earlier 
time t* is just cdt* — R*dx or dx — cdt*/R*. Upon substituting dx into equation (|A6[) we obtain, 



dN c h*(a, z)(R*x) 2 dad[i 
d£l ~ 1 + z 



dt* 
dz 



dz 



(A7) 
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where we have also expressed everything in terms of redshift z by transforming from the variable f* to z. 

An expression for (i?*x) 2 is derived by recalling the relationship between energy flux <&e and luminosity distance 
d,L, namely, 



d£ 
dAdt 



= $K = 



4ird 2 L 4ird 2 L dt* 



d£* , 7? \ — i / d£* 
= ( 4 ^l) -J7T 



4ndl 



d£ 



§e dA 



dt 
dt, 



\dAdt J 



dA 



or 



4nd 2 L 



dA = 1 



(A8) 



where the area element dA of the Robertson- Walker metric is 

dA = (Rx) 2 sm6 

Hence 



dA* 

Ha 



Rx 



dn 



R 2 x 2 dfl 
1 

(1 + z) 2 



(A9) 



where we have utilized equation (IA4p and use <if2* = <if2. Substituting equation (|A9|) into equation (|A8|) for cL4, we 
obtain 



(i?*a;) 2 = 



4 



(1 + *) 



4 " 



For a flat ACDM universe (,Peebleslll998h . 

dz 



dt, 



H (l + z)^/n m (l + z) 3 + fl A , 



(A10) 



(All) 



where fl m and £l\ are the ratios of the energy densities of total mass, dominated by dark matter, and dark energy, 
respectively, compared to the critical density. Substituting equations (|A10[) and (|A11|) into equation (|A7|) gives 



dN c d 2 L (z) h,{a, z) da d[i dz 

dfi = ^0(l + z) 6 V^m(l+^) 3 + "A ' 



(A12) 



for randomly oriented GRB sources. After using the relation that the proper density h*(a,z) increases with redshift 
oc (1 + z) 3 in the comoving density h co (a, z), that is h,{a, z) = (1 + z) 3 h co (a, z), we obtain the result given in equation 
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Fig. 9. — Measured and model redshift distributions for SFR3 (Hopkins & Bcacom 2006). In the left column, from top to bottom, 
we have the pre-Swift cumulative redshift distribution with the assumed range of jet opening angles £>j m i n = 0.05, 6*j max = 7r/2, 0.7, and 
0.4 radians, respectively. The right column has results as in the left column but for Swift data. In each panel, each curve represents the 
cumulative redshift distribution for different jet opening angle power-law indices s, where s = —0.6, —1.2, —1.8, and —2.4, from far left to 
far right, respectively. The assumed gamma-ray energy is = 2 X 10 51 ergs. 
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01234567 
Redshift z 

Fig. 12. — SFR2: Swift cumulative redshift distribution with an assumed range of the jet opening angles Bj m i a = 0.05 and $j max = 
0.7 radians. From top to bottom, the assumed gamma-ray energies are 5 X 10 50 , 2 X 10 51 , and 10 52 ergs, respectively. Each curve represents 
the cumulative redshift distribution for different power-law indices s, where s = —0.6, —1.2, —1.8, and —2.4, from far left to far right, 
respectively. 
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Fig. 13. — SFR4: In the left column from top to bottom, we have the pre-Swift cumulative redshift distribution with the assumed 
gamma-ray energy £* 7 = 5 X 10 50 , 2 X 10 51 , and 10 52 ergs, and an assumed range of the jet opening angles between 6j m i n = 0.05 and 
0jmax = 0.7 radians. The left and right columns show the model distributions for pie-Swift and Swift data, respectively. In each panel, 
each curve represents the cumulative redshift distribution for different indices s, where s = —0.6, —1.2, —1.8, and —2.4, from far left to far 
right, respectively. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 



Opening Angle 9j[rad] 

Fig. 14. — SFR4: pve-Swift cumulative opening-angle distribution with an = 10 52 ergs, and an assumed range of the jet opening 
angles between 0j m i n = 0.05 and 0j max = 0.7 radians. Each curve represents the cumulative opening-angle distribution for different indices 
s, where s = —0.6, —1.2, —1.8, and —2.4, from far right to far left, respectively. 
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Fig. 15. — SFR5: In the left column, from top to bottom, we have the pre-Swift cumulative redshift distribution with the jet opening 
angles ranging from £>j m i n = 0.05 to 0j max = 7r/2, 0.7, and 0.4 radians, respectively. In the right column we have the same thing as in the 
left column but for Swift data. In each panel, each curve represents the cumulative redshift distribution for indices s, where s = —0.6, 
— 1.2, —1.8, and —2.4, from far left to far right, respectively. The assumed gamma-ray energy is £* 7 = 2 X 10 ergs. 
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Fig. 16. — SFR5: In the left column from top to bottom, we have the pie-Swift cumulative opening-angle distribution with the jet 
opening angles ranging from $j m i n = 0.05 to 6>j max = n/2, 0.7, and 0.4 radians, respectively. In the right column we have the same thing 
as in the left column but for the Swift data. In each panel, each curve represents the cumulative redshift distribution for indices s = —0.6, 
— 1.2, —1.8, and —2.4, from far right to far left, respectively. The assumed gamma-ray energy is £* 7 = 2 X 10 51 ergs. 
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Fig. 20. — SFR5: Same as Figure [TBI with £, 7 =4x 10 51 ergs. The solid curves are for s = —1.3. 
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Fig. 22.— SFR6: Same as Figure \W\ with £, T =4x 10 51 ergs. 



Redshift Distribution of GRBs 



27 




0.05 0.1 0.15 0.2 0.25 0.3 0.05 0.1 0.15 0.2 0.25 0.3 
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Fig. 23. — Swift cumulative redshift distribution assuming SFR5 (left column) and SFR6 (right column), with £* 7 = 4 X 10 51 ergs. The 
two bottom panels show the above best-fitted results at z < 1 in more detail. 
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Fig. 24. — Integral size distributions (a) for SFR5 (solid) and SFR6 (dash curve), respectively. The inset is the differential size distributions 
(b) for SFR5 (solid) and SFR6 (dash curve), respectively. The filled circles curve represents the 1024 ms trigger timescale from the 4B 
catalog data, containing 1292 bursts. 
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